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I NTRODUCT I  ON 


The  linear  theory  of  the  propagation  of  a  laser  beam  in  a  random  or 
turbulent  medium,  which  is  valid  for  a  low  power  laser,  has  been  widely 
investigated  [_\ ,22.  However,  a  high  energy  laser  beam  will  modify  the 
medium  by  various  effects,  such  as  inducing  turbulence  in  a  quiescent 
medium  or  changing  the  state  of  previously  existing  turbulence.  Such  a 
modified  medium  will  in  turn  change  the  mode  of  propagation.  Experiments 
have  shown  that  a  laser  beam  can  indeed  induce  turbulence,  and  some 
criteria  for  the  onset  of  turbulence  have  been  developed  C3,4,5_.  In 
the  present  work,  we  shall  develop  a  theory  of  thermal  turbulence  as 
induced  by  a  laser  beam.  Here  the  laser  acts  as  a  source  of  heat  which 
is  deposited  into  the  fluid  medium  by  conduction  and  convection,  due  to 
absorption  and  beam  inhomogenei ty ,  respectively;  in  a  self-consistent 
way,  the  thermal  fluctuations  also  modify  the  heat  source.  The  nonlinear 
equation  describing  the  fluctuations  of  temperature  gives  rise  to  a 
hierarchy  which  needs  to  be  closed.  The  method  of  cascade  decomposition, 
which  was  developed  by  Tchen  to  close  the  hierarchy  and  which  has  been 
successfully  applied  to  hydrodynamic  and  plasma  turbulence,  will  be 
extended  to  laser  induced  turbulence  C6,7H. 

In  the  present  work  we  neglect  the  acoustic  effects,  as  is  permitted  in 
all  problems  of  thermal  oscillations  under  the  Boussinesq  approximation. 
The  buoyancy  effect  which  is  retained  in  the  general  formulation  is 
found  to  be  negligible  in  the  inertia  and  dissipation  subranges  of  the 


spectrum. 


FUNDAMENTAL  EQUATIONS 


The  fundamental  equations  which  describe  the  fluide  motion  in  the  atmos¬ 


phere  in  the  presence  of  laser  heating  are  the  Nav i er-Stokes  equation 
of  momentum, 


=  -Vp  +  V(p  v  V  u)  +  4-  V(p  v  V  u)  -6  g(T-T0) 


(I  ) 


the  equation  of  continuity. 


9  p 

y!p  +  V  (p  u)  =  0 


(2) 


and  the  equation  of  heat, 


2 

+  u VT  =  -  p  V  u  +  k  V  T  +  +  ct  I 


(3) 
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Here  p  is  the  density,  u  is  the  fluid  velocity,  T  is  the  temperature, 
cv  is  the  specific  heat  at  constant  volume,  $  is  the  Rayleigh  dissipation 
function,  I  is  the  intensity  of  the  laser  beam,  and  a  is  the  absorption 
coefficient.  The  buoyancy  force  is  represented  in  Eq.  (I)  by 

(p  -  po)  g  =  -  B(T-To)  g  (4) 

where  g  is  the  acceleration  of  gravity,  and  6  is  the  coefficient 
of  expansion. 

It  is  seen  that  the  effect  of  the  laser  resides  in  (a)  the  heat 
disposition  term  al  in  Eq.  (3)  which  modified  the  temperature 
distribution,  and  (b)  the  buoyancy  force  in  the  momentum  equation  (I), 
modifying  the  velocity  distribution. 

In  problems  of  thermal  convection  and  diffusion,  it  is  customary  to 
introduce  the  Boussinesq  approximation,  which  neglects  the  acoustic  or 
compressibility  effect  in  all  terms  except  the  buoyancy  term  of  Eq.  (I). 
This  approximation  reduces  Eq.  (2)  to 

V-u  =  0  (5) 

We  shall  also  neglect  the  Rayleigh  dissipation  function  in  Eq.  (3).  The 
system  of  Eqs.  (I),  (3),  and  (5)  determine  the  modifications  of  the 
atmospheric  motions  by  the  heat  deposition  from  the  laser.  It  i s  to  be 
remarked  that  in  the  absence  of  such  a  heat  deposition,  the  above  system 
of  equations  degenerates  to  the  one  governing  the  classical  problem  of 
therma I  i nstabi I i ty . 

In  the  present  study  it  is  necessary  to  specify  the  heat  deposition  term 
al  of  Eq.  (3).  We  note  that  I  is  related  to  the  index  of  refraction  n 
on  the  basis  of  photon  conservation,  and  subsequently  is  related  to  the 
temperature  by  an  equation  of  state.  In  principle,  these  relations  can 
be  derived  from  time-dependent  microscopic  models.  However,  their 
governing  time  scale  is  much  smaller  than  the  time  scale  of  the  turbulent 
motion,  and  consequently  a  quas i -stat ionary  approximation,  known  as  the 
adiabatic  approximation,  can  be  used. 

This  amounts  to  assuming  that  the  heat  deposition  is  spontaneous  and 
localized.  This  requires  that  the  approach  to  equilibrium  of  the  heating 
process  is  faster  than  that  of  the  turbulent  process.  This  condition  is 
usually  fulfilled.  Under  such  a  circumstance  we  wish  to  find  a  relation 
analogous  to  that  between  the  energy  flux  and  the  gradient  from  the 
thermodynamics  of  irreversible  processes.  A  satisfactory  derivation  of 
such  a  relation  requires  the  analysis  of  the  equation  of  wave  propagation 
in  a  medium  of  variable  index  of  refraction  and  the  study  of  the  absorp¬ 
tion  or  heat  deposition  during  the  propagation.  A  diffusion  theory  of 
absorption  of  the  laser  beam  intensity  on  the  basis  of  a  stochastic 
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theory  will  be  developed  in  a  later  report.  In  this  paper  we  shall 
present  a  heuristic  argument.  To  that  end,  we  know  that  a  variation  of 
temperature  gradient  6(VT),  which  corresponds  to  a  variation  in 
refractive  index  gradient  6(^h),  will  cause  a  thermal  distortion  of  the 
beam  intensity,  expressed  in  the  form 

6 1  o  V6n 

“f  =  %  D-—  (6a) 


This  is  analogous  to  the  relation  between  flux  and  gradient  in  the 
thermodynamics  of  irreversible  processes,  as  mentioned  earlier.  Here  % 
is  a  scale  length  characteristic  of  the  propagation  pathlength  in  the 
turbulent  medium  (it  should  not  be  confused  with  the  inner  and  outer 
scales  of  turbulence),  and  D  is  a  differential  operator 

D  =  V  +  ttn  I  (6b) 

Formula  (6a)  is  a  variant  of  the  equation  developed  by  Gebhardt  and  Smith 
L8D.  We  may  consider  that  Eq.  (6a)  serves  as  a  relation  of  state  which, 
when  combined  with  the  const i tu i t i ve  equation 


6n 


(7) 


will  give  the  equivalent  value  of  the  temperature  increase  6T  correspond¬ 
ing  to  the  energy  deposition  61,  as  follows: 

61  =  A  D •  V  6  T  (8a) 

with 


We  shall  assume  that  a  turbulent  fluctuation  follows  such  a  variation, 
and  shall  rewrite  Eq.  (8a)  as 

I'  =  A  D*V  6  (9) 

Here  I1  and  6  are  turbulent  fluctuations  of  radiation  intensity  and 
temperature,  respectively,  while  Io  and  To  are  the  corresponding  mean 
quantities;  further. 
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D  =  V  +  V  Jin  I 


(10) 


lo 


1+  i s  +o  be  noted  that  A  and  A  >  0,  because  (9n/9T)p  and  (9n/9T)p  are 
negative  for  most  fluids.  ° 

The  relationship  between  I1  and  0,  as  expressed  by  Eq.  (9),  is  analogous 
to  the  formula  for  the  beam  intensity  developed  by  Gebhardt  and  Smith 
for  the  geometric  optics  regime  [8U.  This  formula  has  been  applied  by 
Collett,  et.  a  I . ,  to  the  determination  of  intensity  correlations  of  a 
laser  beam  in  a  turbulent  medium  The  results  obtained  by  this  method 

are  identical  to  the  more  rigorous  method  of  wave  propagation  used  by  Ho, 
when  the  effects  of  diffraction  are  omitted  C I 0U .  That  approximation  is 
valid,  because  the  optical  length  scale  is  very  small  compared  to  the 
scale  of  the  turbulence. 

The  system  of  Eqs.  (3)  and  (9)  forms  the  fundamental  differential 
equations  governing  the  interaction  of  a  laser  beam  with  the  atmospheric 
medium,  as  described  by  the  variables  u,  p,  T,  and  I1. 

EQUATIONS  OF  TURBULENT  MOTIONS 

We  decompose  the  four  variables  p,  T,  u,  and  I  into  a  mean  part  and  a 
f I uctuation: 

u  =  U  +  u 
p  =  Po  +  p' 

t  =  To  +  e 

I  =  Io  +  I'  (II) 

In  the  absence  of  a  mean  wind,  u  represents  a  velocity  fluctuation. 

By  substituting  Eq.  (II)  into  Eqs.  (I),  (3),  and  (5),  we  find  the  follow¬ 
ing  four  equations  determining  the  four  variables  u',  p',  0,  and  I': 


9U 

3+  +  u  V  u 


1  _ 

"  V  p1  +  v  ^u-0g0  +  u  Vu 


(12) 


V  u  =  0  (13) 

||-  +  u  TO  =  -u  VT0  +  u  W  +  (k/pcv)  V20  +  (a/pcv)  I'  (14) 
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I1  =  A  (  V20  -  r  SB) 


(15) 


with 


r  =  -  V  £n  I0 


(16a) 


where  r-1  represents  the  scale  of  the  laser  beam. 

As  a  result  of  the  turbulent  heating  of  the  atmospheric  medium  by  a  beam, 
we  expect  a  distribution  of  the  intensity  of  temperature  fluctuation, 
with 


A  =  -  V  £n  02  (16b) 


where  A_1  is  a  scale  in  the  profile  of  02  analogous  to  Eq .  (16a). 

1+  is  well  known  that  the  buoyancy  controls  the  large  scale  motions  of 
the  velocity  spectrum.  This  effect  is  negligible  for  the  inertia  and 
dissipation  subranges  under  the  present  investigation,  enabling  us  to 
decouple  Eqs.  (12)  and  (13)  from  Eq .  (14)  and  (15).  Hence,  we  need  only 
to  be  concerned  with  the  thermal  turbulence,  and  regard  the  velocity 
spectrum  as  a  known  quantity.  In  such  a  circumstance,  and  with  a 
substitution  of  Eq .  (15)  into  Eq.  (14),  we  find  the  thermal  turbulence 
governed  by  the  following  equation  of  heat  transfer: 


30  —  9 

3-f-  +  g*V0  -  -  u  •  VTQ  +  V*u0  +  (X  +  y)  V*-0  -  yT»V0 


(17) 


The  left  hand  side  of  Eq.  (17)  represents  the  time  rate  of  change  and 
the  convection  of  the  temperature  fluctuation.  The  right  hand  side 
contains  a  generation  of  temperature  fluctuation  from  the  velocity 
fluctuations  in  the  presence  of  a  mean  temperature  gradient,  a  conduction 
from  the  thermal  flux  by  eddy  motions,  a  molecular  conduction  with  a 
molecular  conductivity 

X  =  k/pcv  (18) 

a  conductivity 

y  =  aA/pcv  (19) 
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from  the  absorption  of  radiation,  and  finally  a  convection  associated 
with  the  beam  in homogeneity. 

CASCADE  DECOMPOSITION 

For  the  determination  of  the  turbulent  spectrum,  it  is  necessary  to  make 
a  Fourier  analysis  of  Eq.  (19).  Instead  of  studying  the  behavior  of 
each  Fourier  component,  we  separate  the  fluctuations  0  and  u  into  two 
groups  and  write 


0 (t,x)  =  0(o) ( t , x )  t  0 ( n (t,x) 


(20) 


and 


u(t,x)  =  u^(t,x)  t  u^(t,x) 


(21  ) 


wheTe  we  define 

,(0) 


ife'  x 


6  (t,x)  =  j dfe'  0(t,fe' )  e  -  - 


and 


0( 1 } ( t , x )  =  dfe'  0 (t ,fe* 


S i mi  I ar I y, 


k 


(22) 


(23) 


u^®^(t,x)  =  -dfe'  u(t,fe')  e 


ife'  x 


(24) 


and 


u(l J(t,x)  =  dfe'  u(t,fe')  e  ~ 


(25) 


where  0^,  0^,  u^  and  u^  are  truncated  within  the  limits  indicated 
by  Eqs.  (22)  to  (25). 
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In  order  to  separate  0^  and  0^  from  0  in  Eq.  (20),  and  u^'  and  ^ 
from  u  in  Eq.  (21),  we  make  use  of  an  ensemble  average 


<■  ■  ■> 


(26) 


which  gives 


<0>(l)  =  e(0),  <e(0)>(l)  =  0<O)f  <e(l))(D  =  0 


(27) 


and 


(28) 


The  decomposition  expressed  by  Eqs.  (20)  and  (21)  implies  that  the 
fluctuations  0  and  u  contain  a  zeroth  rank  and  first  rank,  which  represent 
macroscopic  and  random  fluctuations,  respectively,  with  the  macroscopic 
quantity  being  of  larger  scale  than  the  random  quantity.  The  ensemble 
average,  Eq.  (26),  of  first  rank  is  an  average  over  many  realizations  with 
identical  macroscopic  conditions,  as  represented  by  0(0)  and  u^,  but 
with  random  0(1)  and  u(l).  In  this  manner,  the  average  over  the  random 
quantities  vanishes  and  the  same  average  over  macroscopic  quantities 
remains  unaffected  by  Eqs.  (27)  and  (28). 

The  ensemble  average  of  zeroth  rank 


(29) 


<•  •  *> 


is  equivalent  to  the  global  average  denoted  by  a  bar  in  a  homogeneous 


turbu I ence. 

By  means  of  the  cascade  ensemble  average,  Eq.  (26),  we  can  screen  out  the 
ranks  0(0)  and  q(I)  from  0  in  Eq .  (20),  and  hence  transform  the  equation 
of  evolution  of  0,  Eq .  (19),  into  correspond i ng  equations  of  evolutions 
of  0(0)  and  0(1).  The  details  for  this  type  of  calculation  can  be  found 
in  the  work  of  Tchen  and  will  be  omitted  here  The  result  of  this 

ca Iculation  leads  to 


u(0)-vto  +  (x  +  p)v20(o)  -  pr-v0(o)  -<\j( 1 V0(  1  *>( 1 ) 


(30) 


and 


(31  ) 


=  -u(  1  ^.V(To+0(O)  )  +  (A+y  )V20(  1  )-yT-V0(  1  J+<u(  1  V0(  1  1  } 

noting  that  the  addition  of  Eqs.  (30)  and  (31)  recovers  Eq.  (19).  The 
last  term  on  the  right  hand  side  of  Eq.  (31)  will  be  omitted  in  the 
following,  since  it  does  not  contribute  to  the  formulation  of 
corre I ati ons. 

EQUATION  OF  SPECTRAL  BALANCE 

We  multiply  Eq.  (30)  by  0'®^  and  take  an  ensemble  average  of  zeroth  rai 
as  defined  by  Eq.  (29).  We  then  find 


=  -<u(0).e(0)>(0)  VT^  _  ( 0 ) v .  ^ ( I  )  ed  )>d  )>(0) 


(O).fl(OK(O) 


+  (A+y )  ^0 


(0) v20 (0) >(0) 


(32) 


or,  rewritten  as 


-(A+y)  <(V0(O))2VU)  +  {j(Aty)V2<0 


(0) \2\(0) 


(0)2v(0)  j 


(33) 


where 


A(0)  E  V  £n  <^e(0)2v(0) 


(34) 
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In  order  to  simplify  Eq.  (33),  we  assume  that  the  turbulent  motions  are 
locally  homogeneous.  This  implies  that  the  spectral  distributions  are 
similarly  dependent  on  position  through  the  intermediary  of  the  para¬ 
meters,  which  are  given  by  the  mean  quantities,  e.g.,  profiles  of  beam, 
mean  temperature  and  intensity  of  thermal  turbulence,  and  the  rate  of 
dissipation  of  the  thermal  turbulence.  Under  that  hypothesis,  we  can 
neglect  the  terms  between  {...}  in  Eq.  (33).  The  term  proportional  to 
T-A(O)  is  a  coupling  between  the  laser  beam  and  thermal  turbulence,  and 
is  therefore  proportional  to  the  product  of  the  gradients  in  the  beam 
profile  and  the  turbulent  intensity  profile.  In  this  connection,  we 
approximate  A^O)  -  A,  where  A  is  defined  by  Eq.  (16b),  since  the  rank  of 
Eq.  (22)  up  to  a  wavenumber  k  lying  in  the  inertia  subrange  embodies  the 
major  portion  of  the  thermal  energy. 

With  the  above  approximations,  we  reduce  Eq.  (34)  to  the  following  form: 

<e(t»2><0)  = .  Cu<o>e«»><o>n=  yiyiyoyo) 

-(X+p)  <(V0(O)  )2>(0)  +  y<0(O)2>(O)  (35a 


with 


Y  =  p  T*A  (35b) 

We  conclude  that  the  rate  of  change  of  thermal  energy  <0  is 

governed  bv  a  production  function  -<u^®  0  ®  VTo,  a  transfer  function 

- (u d  ^ 0 ^  ^0 (0)  y  0) ,  a  dissipation  function  -(X+p ^(VO (0)  )2^(0)  t 

combining  the  effects  of  molecular  conduction  and  the  absorption  of 
radiation,  and  finally  a  source  function  y<9^)2>(0)  from  the  laser  beam. 

HEAT  FLUX  AND  EDDY  CONDUCTIVITY 

In  Eq.  (35)  there  occurs  a  term, 

<y( 1 )0( 1 )>( I )  (36) 

called  eddy  "thermal  flux"  of  first  rank,  and  a  similar  term  of  zeroth 
rank. 


<u(O)0(O)>(O)  =  <U0>  -  <u(l >0<l)>(l ) 


(37) 
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The  statistical  study  of  molecular  transport  processes  can  be  based  upon 
a  Langevin  equation,  which  describes  the  time  evolution  of  a  fluctuation 
in  the  Lagrangian  representation.  We  shall  generalize  the  above  Langevin 
equation  to  apply  to  a  turbulent  transport.  In  order  to  do  this,  we  note 
that  Eq.  (31)  may  serve  this  purpose  by  considering  the  variable  t  in  the 
Lagrangian  time  derivative  as  a  single  variable,  since  this  is  the 
requirement  for  the  calculation  of  the  transport  properties.  Effectively, 
we  regard  the  Fourier  form  of  Eq.  (31)  as  a  turbulent  Langevin  equation 
with  k  as  a  parameter.  It  will  be  written  as  follows: 


d6(l  }(t,fe)  =  _  ,,(|  ) 

dt 


(t,fe)  v(Tote(0)) 


x*fe2e ( 1 > (t,fe) 


(38) 


where  T0  and  0^  have  weaker  inhomogeneities  than  6^  and  u^,  and 
\*k 2  represents 


\*k2  =  Xk2  +  u(fc2  +  r-A) 


(39) 


A  formal  sol ution  of  Eq .  (38)  is 

t 

e(  1  > (t,fc)  =  -  V(To+e(0) )  •  dt'  expC-A*k2 Ct-t' )D  i/^(t’,fe) 

j 


+  0^^(t=o,fe)  exp(-A*fe2t) 


(40) 


Here  we  have  put  V(T0  +  0^))  outside  of  the  integral,  because  it  varies 
slowly  with  time.  It  follows  from  Eq.  (40)  that 


<uU)(t,fe)  1  ‘  (t,fe))  =  -  nv'.  #(fc) 

J 


,(D 


(  I  ) 


( I  )  ,  3 (To  +  0(O)) 


3xc 


(41  ) 


where 


(I  ) 


I;(fe)  =  :  dt  ’  expC~A*fe2  (t-t' )H  /u^^(t',k)  uV;(t,-fe)\ 

Sj  '  J  \  s  ~  j  -  / 

O 

noting  that 

<0(l  Uo,k )  0(l  )(t,-fe))(l  }  =  0 


(I  ) 


(I) 


(42) 


(43) 
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for  t  larger  than  the  correlation  time.  By  the  same  token  we  can  replace 
the  upper  limit  t  by  ~  i n  Eq.  (42)  to  give 


nU)(fe>  - 

sj  ' 


dt’  exp  EA*fe2(t-t’)J  <u^(t',fe)  uU)(t,  -k)') 


(44) 


O 


For  an  isotropic  eddy  conductivity  we  have 


(45) 


and 


I  ) 
n 


Uz’t 


r 

^  df  expCA*fe2(t-t')D  <u(l)(t’,fe)  u(l)(t,-fe)>(l)  (46) 


This  reduces  the  thermal  flux,  Eq .  (36),  to 
<u!n(t,-fe)  e(l  }(t,fe)>(l  ]  =  -  n(l)(fe) 


3(T0  +  9(0)) 

3XJ 


(47) 


i n  k  space,  or  to 


<uU)(t,x)  eu)(t,x)>(l } 
J 


-  n 


( i )  3 (To  +  e(0)) 
3xj 


(48) 


where  n^ * ^  is  the  eddy  conductivity  in  x  space. 

The  turbulent  transport  relations,  Eqs.  (41)  and  (48),  take  the  form 
analogous  to  the  phenomenological  relations  of  Onsager  in  the  thermo¬ 
dynamics  of  irreversible  phenomena,  stating  that  any  flux  is  caused  by 
the  contribution  of  a  force  or  gradient  with  a  proportionality  constant 
called  the  transport  coefficient.  Here  we  have  derived  such  relations  by 
means  of  a  generalized  Langevin  equation,  Eq.  (38),  applicable  to 
turbulent  motions;  hence,  we  have  determined  the  relevant  forces  and 
transport  coefficient,  as  given  by  Eqs.  (41).  (44),  (46),  and  (48).  A 
precise  determination  of  the  structure  of  t/  U  jn  terms  of  the  spectral 
function  F(k),  such  that 

k 

^  ^(0)2^0)  =  dfet  F(fet  )  (49) 
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is  given  by  Tchen  C6D*  We  shall  not  enter  into  its  detailed  derivation, 
but  simply  write  the  following  result: 


n' 


I  dfe» 

k 


F(fe'  ) 
w^  ^  (fe1 ) 


Here  the  relaxation 

w(D(fe')  =  fe«2 

=  fe'2 


frequency  is 
(X*  +  1 ^ ) 


(50) 


(51 ) 


since  the  molecular  conductivity  X*  can  be  assumed  to  be  smaller  than  the 
eddy  conductivity.  The  integral  equation,  Eq .  (50),  with  the  relaxation 
frequency  given  by  Eq.  (51),  leads  to  the  following  solution: 


.(I)  _ 


2  c. 


i dk'  F(fe')  fe'  2 


T  1 

2 


(52) 


where  C2  is  a  numerical  coefficient  evaluated  for  three-dimensional 
turbulence  to  be  Q6D 


c2 


(53) 


EQUILIBRIUM  RANGE  OF  THE  SPECTRUM 


Along  with  the  spectral  fluctuation  distribution  F(fc)  for  the  velocity, 
Eq.  (49),  we  introduce  a  spectral  distribution  G(fc)  for  the  temperature 
which  is  defined  to  be 


^<e<0)2>  <°>  = 


k 

r 

dfe'  G(fe') 


(54) 


The  use  of  Eq.  (48)  permits  the  transfer  function  in  Eq .  (35)  to  be 
rewritten  as 
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<v-<u(  1  V  1  )x/ 1 )  e(0)>(0)  =  _  n(i)  <V0)v2(to  +  e(0)}>(0) 

=  _  n(  i )  ^e(0)v2e(°)'/(0) 

=  n(  1 J  ((.  ve(0)  )2>e(0)-  1  n(  1  )v2<^e(0)2>(0) 

=  n(  1  ]  <(Ve(0))2;,(0)  (55) 

with  the  application  of  local  homogeneity  of  turbulence. 

If  we  introduce  the  notations 

fe 

jo  =  <^(V0(O))2  \(0)  =  2  dfe'  fe'2  G(fe')  (56) 


o 

and 


(0)  .  ( I  ) 

n  =  n  +  n 

we  can  rewrite  the  equation  of  spectral  balance,  Eq..(35a),  in  the  form: 

fe 

jp  dfe'  G(fe')  =  -  (n-nU))  (VT0)2  -  Q+y)j°  -  n(l  J°  (58) 

/ 

°  fe 
+  y  dfe'  G ( fe '  ) 


or 


■|p  dfe'  G(fe")  =  (n  -  n(  1 )  )J 
°  fe 

+  y  dfe'  G(fe' ) 

O 


-  Q0  +  n  )J 


(59) 


where 
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J  =  (VTo)2 


(60a) 


and 


XQ  =  X  +  p 


C60b ) 


are  called  thermal  vor+ici+y  and  total  conductivity,  respectively.  The 
total  conductivity  includes  a  molecular  conductivity  X  and  a  conductivity 
p  enhanced  by  the  beam.  Use  has  been  made  of  Eqs.  (56)  and  (57). 

For  k  =  °°,  we  reduce  Eq.  (59)  to 


at 


dfe!  GCk1 )  =  nJ  “ 


XQ  J 


+  Y  dk'  G(k' ) 


o 


o 


(61)) 


as 


n(l )(fe=-)  =0,  J  =  J°(k=~)  (62) 

The  equilibrium  range  of  the  spectrum  is  defined  as  that  range  covering 
sufficiently  large  wave  numbers  such  that,  in  the  integral  term  on  the 
left  hand  side  of  Eq.  (59),  the  upper  limit  k  may  be  replaced  by  °°. 

Under  that  circumstance,  and  by  subtracting  Eq .  (59)  from  Eq.  (61),  we 
obtain  the  equation  of  spectal  balance  in  the  equilibrium  range, 

CO 

n(l)T  +  (Xo+n(l))J°  +  y  dk'  G(k')  =  XoJ  (63) 

k 

THE  SPECTRUM  IN  THE  ABSENCE  OF  BUOYANCY 

The  universal  range  of  the  spectrum  is  usually  divided  into  the  following 
subranges  in  the  increasing  orders  of  wavenumbers:  buoyancy,  inertia, 
and  dissipation.  As  mentioned  earlier,  the  buoyancy  subrange  reigns  over 
the  largest  scale  of  the  spectrum,  and  will  not  be  included  in  the 
present  study.  Consequently,  the  velocity  spectrum  in  the  inertial  sub¬ 
range  is  given  as  the  Kolmogorov  law, 
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F  =  A  e2/ 3  r5/3 


C64a ) 


Here  e  is  the  rate  of  energy  dissipation  and  the  numerical  coefficient 
A  i s  C63 


A  =  1.58 


(64b) 


The  spectrum,  Eq.  (64a),  is  the  most  commonly  adopted  law  for  the 
description  of  the  interaction  of  a  propagating  laser  beam  in  a  turbulent 
atmosphere.  Following  that  law,  we  find  the  eddy  conductivity  from  Eq. 


(52)  to  be 


n(l)  =  c3  e1/3  rV3 
with  C6] 


(65a) 


,  I 

c„  =  (_  Ac-?) 2  5  1.26 
3  4  z 


(65b) 


For  the  derivation  of  the  spectral  function  G(fe)  for  the  temperature 
fluctuations,  we  consider  Eq.  (63),  which  we  rewrite  in  the  form 


(Ao  +  n( 1 } )  (j  +  j°)  +  y 


Jdfe' 


G (fe ’ )  =  Ao ( J  +  J  ) 


or,  upon  dividing  by  A0  +  n^. 


(66) 


J  +  J °  +  "y ( A o  +  n(  1 } ) 


dfe'  G(fe' )  =  A0 (Ao+n( 1  5 


l  _ 

(J  +  J) 


fe 


(67) 


We  shall  introduce  the  following  notations: 

oo 

H(fe)  =  !  dfe'  G(fe' )  (68a) 

k 

P  =  -r^-4r  In  (X o  +  n(l  5)  (68b) 

l-$  dfe 
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J 


A°  ( J  +  J) 


i.n  (Ao  +  n(  1  * ) 


(68c) 


0  = 


2  ( Ao  + 


n(l >)fe2 (I-*) 


d_ 

dfe 


$  =  - 1 - 

2fe2 ( Ao  +  n(l >) 


(68d) 


and  transform  the  integral  equation,  Eq.  (67),  into  the  following 
differential  equation 


dH 

dF 


+  P 


H  =  Q 


(69) 


with  the  solution 


G(fe )  =  dfe'  Q(fe')  exp[-  dfe"  P(fe")D 


=  -  Q(fe)  -  P(fe)  dfe'  Q(fe')  exp[-  dfe"  P(fe")] 


or,  in  an  approximate  form. 


(70a) 


G(fe)  =  -  Q(fe)  -  P(fe)  dfe'  Q(fe')  (70b) 

fe 

valid  for  a  weak  beam  factor  y,  or  4>. 

Equation  (70b),  which  is  valid  for  all  subranges,  appears  in  a  form  too 
complex  for  interpretation.  Therefore  we  shall  consider  them  separately. 
We  find: 

(a)  In  the  inertia  subrange 

G(fe)  =  J-  e,  e”1/3(l-*.)  fe"5/3  (71) 

3c  -j  A  1 
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with 


4>j  =  (2c-j)  '(y3/e)1//3  fe-2^3  (72) 

(b)  In  the  dissipation  subrange 

G(fe)  =  el/3  e^X“2  (  l+<J>d  )k  ~  13/3  (73) 

with 

<t >  =  (y/ 9Xo )k~2  (74) 

d  z 

In  the  above,  e  and  are  rates  of  dissipations 


e  =  v  (Vy)2  (75) 


ex  =  X0  (V9)2  =  XQJ 


(76) 


for  the  fluid  and  temperature  fluctuations,  respectively. 

The  transition  from  the  inertia  to  the  dissipation  subranges  occurs  at  a 
critical  wavenumber 

kX  =  c33/4(£/*o)l/4  (77 ) 


We  conclude  that  the  effect  of  the  beam  factor  is  to  decrease  the  power 
of  -5/3  in  the  inertia  subrange  and  steepen  the  power  of  -13/3  in  the 
dissipation  subrange.  The  laser  heating  increases  Xo  and  J. 

CONCLUSIONS 

In  the  present  work  we  have  shown  that  a  laser  can  induce  temperature 
fluctuations  which,  because  of  their  nonlinear  nature,  can  lead  to 
turbulence.  Because  the  hydrodynamic  system  is  nonlinear,  it  generates 
a  "hierarchy,"  i.e.,  an  infinite  sequence  of  equations  of  ever  increasing 
order  of  correlations.  This  is  closed  by  the  method  of  cascade. 
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The  universal  subrange,  which  is  divided  into  an  inertia  subrange  and 
a  dissipation  subrange,  is  not  controlled  by  the  buoyancy  effect,  so 
that  the  spectrum  of  velocity  fluctuations  follows  the  Kolmogorov  law. 

The  spectrum  of  temperature  fluctuations  is  derived  for  the  inertia  and 
dissipation  subranges,  and  includes  the  effects  of  enhanced  diffusion 
and  inhomogeneity  of  the  laser  beam.  The  resulting  equation  of  spectral 
balance  is  in  the  form  of  a  nonlinear  i ntegro-d i f ferent i a  I  equation  and 
determines  the  spectral  distribution  of  the  temperature  fluctuations; 
it  is  solved  for  a  turbulent  state  in  statistical  equilibrium.  It 
contains  a  production  function  due  to  the  beam  i n homogene i ty,  a  transfer 
function  which  describes  the  transfer  of  energy  across  the  spectrum 
from  large  eddies  into  small  eddies,  and  finally  a  dissipation  function 
which  determines  the  dissipation  of  small  eddies  by  molecular  motions. 

If  the  temperature  variable  is  considered  as  a  passive  scaler  driven  by 
a  fluid  turbulence  obeying  a  Kolmogorov,  it  is  known,  by  dimensional 
theories,  that  the  temperature  fluctuation  will  also  follow  the 
Kol mogorov | aw  in  the  inertia  subrange,  and  a  fe-13/3  law  in  the  dissipa¬ 
tion  subrange.  By  means  of  a  cascade  theory,  we  have  derived  analyti¬ 
cally  the  temperature  spectra  which  confirm  the  results  of  the  dimensional 
theories.  The  heating  by  a  laser  beam  introduces  the  following  two 
effects : 

(1)  In  view  of  the  increased  conductivity  by  the  absorption,  the 
rate  of  thermal  dissipation  £\  increases,  and,  consequently,  the 
temperature  spectrum  has  a  higher  value  in  the  inertia  subrange,  and  a 
lower  value  in  the  dissipation  subrange.  For  the  same  reason,  the 
critical  wavenumber  separating  the  two  subranges  occurs  at  a  lower  value. 

(2)  The  laser  beam  has  an  inhomogeneous  profile,  denoted  by  the 
beam  factor  y  (see  Eq.  (35b)).  Since  the  inhomogeneity  is  usually  weak, 
the  powers  -5/3  and  -13/3  are  modified  slightly,  with  a  negative  and  a 
positive  correction  term  in  the  inertia  and  dissipation  laws, 
respectively  (see  Eqs.  (71)  and  (73)). 
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ATTN:  CODE  2627 

Washington,  D.  C.  20390 
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Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-XL-D 

Ft  Monmouth,  NJ  07703 

CO,  Edgewood  Arsenal 

Defense  Devel  <Sr  Engineering  Lab 

ATTN:  SMUEA-DDW  2,  Mr.  H  Tannenbaum 
Edgewood  Arsenal,  MD  21010 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-WL-D 

Ft  Monmouth,  NJ  07703 

Commanding  Officer 

Aberdeen  Proving  Ground 

ATTN:  STEAP-TL 

Aberdeen  Proving  Ground,  MD  21002 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-NL-D 

Ft  Monmouth,  NJ  07703 

Director,  Ballistic  Rsch  Labs 

US  Army  Aberdeen  Rsch  &  Dev  Cen 

ATTN:  AMXRD-BTL,  Mr.  F  J.  Allen 
Aberdeen  Proving  Ground,  MD  21005 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-KLD 

Ft  Monmouth,  NJ  07703 

Commanding  Officer 

US  Army  Ballistic  Rsch  Labs 

ATTN:  AMXBR-B  &  AMXBR-1A 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-VL-D 

Ft  Monmouth,  NJ  07703 
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Commanding  Officer 

US  Army  Land  Warfare  Lab 

ATTN:  CRDLWL-4A 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-CT-D 

Ft  Monmouth,  NJ  07703 
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Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer 

US  Army  Elect  Proving  Ground 

ATTN:  STEEP  MT-EA 

Ft  Huachuca,  Arizona  85613 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-SC 

Ft  Monmouth,  NJ  07703 

Commanding  Officer 

US  Army  Dugway  Proving  Ground  Library 
ATTN:  STEDP-TL,  Tech  Library 

Dugway,  Utah  84022 

Commanding  General 

US  Army  Elect  Command 

ATTN:  AMSEL-RD-D 

Ft  Monmouth,  NJ  07703 

Commanding  Officer 

US  Army  Cold  Regions  R&E  Lab 

ATTN:  CRREL-RP  (Dr.  Yin  Chao  Yen) 
Hanover,  New  Hampshire  03755 

Commanding  General 

US  Army  Combat  Dev  Cmd 

ATTN:  CDCMR-E 

Ft  Belvoir,  VA  22060 

Commanding  Officer 

Yuma  Proving  Ground 

ATTN:  STEYP  AD  (TECH  LIBRARY) 

Yuma,  Arizona  85364 

President 

US  Army  Artillery  Board 

Ft  Sill,  OK  73503 

Commanding  Officer 

US  Army  Arctic  Test  Center 

APO  Seattle  98733 
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Dir  of  Defense  Rsch  &  Engin 

ATTN :  Tech  Library 

RM  3E-1039,  The  Pentagon 

Washington,  D.  C.  20301 

CO,  US  Army  Tropic  Test  Center 

ATTN :  STETC  AD-TL 

APO  New  York  09827 

2 

Defense  Intelligence  Agcy 

ATTN:  DIAAP-10A2 

Washington,  D.  C.  20301 

Commanding  Officer 

USACDC  Inst  of  Nuclear  Studies 

Ft  Bliss,  TX  79916 

Commanding  Officer 
USADDC  Inst  of  Systems  Analysis 
ATTN:  CDISA-SSCI  &  CDISA-STANO 
Ft  Belvoir,  VA  22060 

Commanding  General 
US  Army  Missile  Command 
ATTN:  AMSMI-RFG,  Mr.  N.  Bell 
Redstone  Arsenal,  AL  35809 

Commanding  Officer 
Frankford  Arsenal 

ATTN:  SMUFA-N-5300,  Mr.  M.  Schoenfield 
Bldg  201 

Philadelphia,  PA  19137 

Commanding  Officer 
US  Army  Satellite  Comm  Agcy 
ATTN:  AMCPM-SC-3 
Ft  Monmouth,  NJ  07703 

Director 

Night  Vision  Lab  (USAECOM) 

ATTN:  ASMEL-NV-OR,  Mr.  S.  Segal 
Ft  Belvoir,  VA  22060 

Sylvania  Elect  Sys-Western  Div 
ATTN:  Technical  Reports  Lib 
P.  O.  Box  205 
Mountain  View,  CA  94040 

NASA  HEADQUARTERS 
Prg.  Ch,  MET  &  Soundings  (SRM) 

Earth  Observations  Programs 
Washington,  D.  C.  20546 

Library-R51-Tech  Reports 
Environmental  Research  Labs 
NO  A  A 

Boulder,  Colorado,  80302 


Commandant 

US  Army  Air  Defense  School 
ATTN:  C&S  Dept,  MSL  SCI  DIV 
Ft  Bliss,  TX  79916 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station  (Bldg  5) 

Alexandria,  VA  22314 
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Director,  USN  Research  Lab. 

ATTN:  Code  6530,  H.  Shenker 
Washington,  D.  C.  20390 

Mr.  E.  F.  Corwin 
Head,  Met  Research  Branch 
Meteorological  Division 
Naval  Air  Sys  Cmd  (AIR-5401) 

Washington,  D.  C.  20360 

Commanding  Officer 
US  Army  Combat  Dev  Cmd 
Intelligency  Agency 
ATTN:  CDCINT-P,  Cpt  Thoresen 
Ft.  Huachuca,  Arizona 

Director 

US  Army  Advanced  Materiel  Concepts  Agency 
2461  Eisenhower  Avenue 
Alexandria,  VA  22314 

University  of  Oklahoma 
Research  Institute 
Field  Artillery  Research  Off 
P.  O.  Box  3124 
Ft.  Sill,  OK  73503 

Dir  US  Naval  Research  Lab 
ATTN:  CODE  6530,  H  Shenker 
Washington,  D.  C.  20390 


